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This paper reviews the current work in the field of Schumann resonances and discusses the topics 
of waveforms and frequency estimates, resonance frequencies and Q factors, source distributions and 
noise spectra, diurnal and seasonal changes of power spectra, and variations of resonance frequencies. 
Observations of the above phenomena are correlated with theoretical considerations. 



1. Introduction 

The resonances in the spherical shell between the 
earth and ionosphere that are due to azimuthal waves 
have been considered in two recent survey papers 
[Galejs, 1964a; Wait, 1964a]. This work contains 
experimental data as well as the required theory. 
The basic theory of the resonances has been reviewed 
by Wait [1964b]. In the present summary the 
emphasis will be on correlation between measurements 
and available theory. There will be no extensive 
analytical developments, but frequent references will 
be made to literature. 

Based on theoretical considerations, Schumann 
[1952a, b, 1957] postulated resonances of the earth- 
to-ionosphere cavity. Koenig [1958, 1959, 1961] 
obtained the first experimental indication of Schumann 
resonances by observing noise waveforms in the out- 
put of a narrow band amplifier. Detailed frequency 
spectra of this noise were first obtained by Balser and 
Wagner [I960]. Other measurements have been 
reported by Fournier [I960]; Benoit and Houri [1961, 
1962]; Lokken et al. [1961, 1962]; Polk and Fitchen 
[1962]; Gendrin and Stefant [1962a]; Balser and Wagner 
[1962a, b, 1963]; Rycroft [1963]; Chapman and Jones 
[1964]. 

The fundamental theory of Schumann resonances 
[Schumann 1952a, b, 1957] is discussed in a book by 
Wait [1962]. Raemer [1961a, b] considers the observ- 
able noise spectra as the response of the earth-to- 
ionosphere cavity due to lightning flashes all over the 
world, but the homogeneous sharply bounded iono- 
sphere model of Raemer introduces high losses, and he 
does not succeed in reproducing the spectral measure- 
ments of Balser and Wagner [I960]. Galejs [1961a, b] 
uses an isotropic ionosphere model of exponentially 
increasing conductivity which is based on measured 
or calculated characteristics of the lower ionosphere. 

1 Paper presented at the ULF Symposium, Boulder, Colo., Aug. 17-20, 1964. 



This model permits a close reproduction of the noise 
spectra measurements in the resonance region, and 
it also provides an agreement with measured ELF 
attenuation rates [Chapman and Macario, 1956; 
Jean et al., 1961]. 

Models of an anisotropic ionosphere are difficult to 
apply to the problem of earth-to-ionosphere cavity 
resonances because of the variations of the magnetic 
field vector along the surface of the earth. The prop- 
agation parameters can be estimated in the presence 
of a horizontal magnetic field component from the 
work of Galejs and Row [1964] and of Galejs [1964b]. 
The earth-to-ionosphere cavity resonances have been 
analyzed for a nonhomogeneous ionosphere in the 
presence of a radial (vertical) magnetic field by Thomp- 
son [1963] and Galejs [1965]. A thin-shell approxi- 
mation of the ionosphere with a superimposed radial 
magnetic field has been considered by Wait [1964c]. 
Some other work of Wait [1963a, b], although intended 
for a different frequency range, can also be applied 
to the earth-to-ionosphere resonance problem in the 
presence of a radial magnetic field. 

Additional references are available in the compre- 
hensive bibliography by Brock-Nannestad [1962]. 

2. Waveforms and Frequency Estimates 

In wide band recordings of atmospheric signals an 
ELF component or slow tail follows the initial VLF 
component of the transient [Hepburn and Pierce, 
1953; Liebermann, 1956a, b: Tepley, 1959; Pierce, 
I960]. An oscillatory structure of the signal becomes 
noticeable by passing it through a band-pass filter. 
An example of a noise recording in the output of a 2 to 
30 c/s band-pass filter [Lokken, 1964, private communi- 
cation] is shown in figure 1. The two horizontal mag- 
netic field components have been observed simul- 
taneously in Canada and in Antarctica. There is a 
correlation between the high intensity ELF noise bursts 
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observed at the various sites. A recording of one of 
the waveforms is shown in figure 2. The darkened 
areas indicate resonance frequencies around 8, 14, 
and 20 c/s. Noise waveforms in the 5 to 20 c/s band 
have been reported by Polk and Fitchen [1962] and 
Polk [1962]. Figure 3 shows sample recordings of 
the outputs of two coils the axes of which are oriented 
in the north to south (N-S) and east to west (E-W) 
directions. Frequency of these waveforms has been 
estimated by counting the cycles between the second 
markers shown on the records. The average fre- 
quency computed in this procedure should be higher 
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FIGURE 3. Noise recording in the 5 to 20 c/s band. 
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FIGURE 1. Noise recording in the 2 to 30 c/s band. 
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Figure 2. Recording of the noise waveform X ofWestham Island, 
B.C. 



than the resonance frequencies of a spectral analysis 
[Galejs, 1962a, 1964a]. Noise recordings made simul- 
taneously in Germany and in Rhode Island have been 
reported by Keefe, Polk, and Koenig [1964]. 



3. Resonance Frequencies and Q Factors 
3.1. Definitions 

A resonance frequency and Q factor may be deter- 
mined from spectral measurements. Analytically it 
is also possible to compute the power spectrum of the 
noise waveform. However, it is more convenient to 
relate the resonance frequency and the Q factor to 
the analytically determined propagation parameters. 
The radial electrical field due to the vertical dipole 
excitation can be represented by 



Er = ~ 



ilds v{v-\- 1) 
^coeha 2 sin vtt 



P v (- cos 0) 



/here 



v(v+l) = (k aS) 2 . 



(1) 



(2) 



For a wave that propagates in the spherical shell be- 
tween the earth and the ionosphere S is defined as the 
ratio of the complex wave number k to free space 
wave number k (S = kjko). The real part of S can be 
seen to be inversely proportional to the phase velocity, 
and the imaginary part of S is proportional to the at- 
tenuation constant. This parameter S is determined 
from solving the appropriate boundary value problem 
of the earth, the air space and the chosen ionosphere 
model [Wait, 1962; Galejs, 1964a]. The resonance 
frequencies of the nth mode are determined from the 
minima of sin vtt in (1) [Galejs, 1962a, 1964a] which 
gives 



fn 



7.5Vra(re+l) 
ReS 



(3) 



Typical values of ReS are 1.4 to 1.2 for frequencies in 
the resonance region. 
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The cavity Q may be determined as a ratio between 
the stored energy and energy loss per cycle or simply 
by the width of the resonance curves. As shown in 
the appendix, in the present problem the electrically 
stored energy W e differs from the magnetically stored 
energy W h , and the Q definition should consider the 
sum of the average stored energies W<> and W tl . This 
expression for Q is 



Q = 



(ReS) 2 + 1 
4ReS ImS' 



(4) 



When considering only the magnetically stored energy 
Q is given by 



Q= 



1 



2 ReS ImS 



(5) 



For a given set of propagation parameters this will 
give a lower Q figure than (4). This expression has 
been used in the past work of Galejs [1961a, 1962a, b, 
1964a]. Considering only the electrically stored 
energy gives 



Q = 



ReS 
2lmS' 



(6) 



The last expression can also be derived from the con- 
cept of complex resonance frequencies [Wait, 1964a] 
or by considering the half-power band width of the 
resonance curve. All three Q definitions yield the 
same results as ReS approaches unity, but they will 
differ near the lower resonance frequency of the earth 
ionosphere cavity. Equation (4) will be used in the 
present calculation. 

Experimentally determined Q values that are esti- 
mated from the half-power bandwidths of the reso- 
nance curves neglect the effects of adjacent resonances 
and of near-field noise that will tend to add to the 
background level. Hence, the half-power level is 
estimated too low, the apparent half-power bandwidth 
is larger and the estimated Q factor may be too low. 

3.2. Characteristics of the Boundary 

In the theoretical models that are advanced for ex- 
plaining the resonance frequencies and the ^-factor 
measurements [Wait, 1962; Galejs, 1964a; Wait, 1964a] 
it is necessary to consider the boundary properties 
of the spherical shell. 

In the frequency range of the resonances the dis- 
placement currents of ground are negligible and the 
ground conductivity cr y =10~ 3 to 1 mho/m is much 
higher than the effective ionospheric conductivity 
o-j that is estimated to be of the order of 10 -7 to 10~ 6 
mho/m. The ground conductivity cr g may be assumed 
as infinite in the first approximation and only the iono- 
spheric properties require a detailed analysis. 

The ionosphere acts as an anisotropic conductor 
and for a z-directed magnetic field its conductivity 



tensor is of the form 



w- 



CTi 







(Jl CT\ 

(To 



(7) 



Equations for computing the conductivity components 
are available in literature [Galejs, 1964a and b] and 
typical conductivity profiles are shown in figure 4 but 
some uncertainties still apply to such models [Yan- 
Zandt, 1964] . The above computations consider the 
ion effects, but assume the operating frequency to be 
zero. These profiles can be further approximated by 
straight lines through the lower ionosphere in the 
semilogarithmic plots of figure 4, which corresponds to 
an expoential height variation of the conductivity com- 
ponents. In a further simplification the anisotropy 
may be ignored (o-i = o- , o- 2 = 0). Even a simple one 
or two layer model can be assumed, such as shown in 
figure 5. These particular models have been used in 
past work by Raemer [1961a and b], Galejs [1962a 
and b], and by Chapman and Jones [1964]. In the 
isotropic exponential model the conductivity is given 
by 



o-(z)=cr(zo) exp [fl(z — z )] 



(8) 



where z = 60 km and the day and night models are 
characterized by cr(zo) =4.63 X 10 -8 mho/m and 
13 = 0.308 km" 1 and by cr(z () ) =6.5 X 10~ 1() mho/m 
and /3 = 0.44 km -1 , respectively. 
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Figure 4. Components of the tensor conductivity F = 0. 
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FIGURE 5. Models of ionospheric conductivity. 



3.3. Single Layer Ionosphere 

The single layer ionosphere model of Raemer [1961a 
and b] assumes the boundary height of the ionosphere 
to be constant, but the conductivity to change with 
frequency so as to obtain the correct resonant fre- 
quencies. However, the resulting Q factors are 
approximately two times lower than in measurements. 
The same difficulties are experienced with an iono- 
sphere model of constant conductivity, but frequency 
dependent boundary height. Similar results are 
obtained by Madden [1961] with a conductivity of 
cr=1.8XlO~ 6 mho/m and the ionospheric boundary 
at h = 80 km. The deficiencies of a single layer 
ionosphere model are less pronounced if it is attempted 
to achieve only a partial match of the resonant fre- 
quencies, as seen from Wait [1964a]. However, 
the Q factors will be too low in particular at the higher 
resonance frequencies (/— 30 c/s). 

Compressibility effects of the ionosphere have been 
considered by Seshadri [1965b] using the model of a 
perfectly conducting flat earth and a sharply bounded 
homogeneous isotropic and lossless ionosphere. The 
compressibility of the ionosphere does not affect the 
wave propagation below it for frequencies in the range 
of Schumann resonances and in the ELF band. How- 
ever, the compressibility should be considered for 
waves that propagate inside the ionosphere [Seshadri, 
1965a]. 



3.4. Two-Layer Ionosphere Models 

Two-layer models give more accurate values of the 
resonance frequencies and of cavity Q or of the com- 
plex propagation parameter S. Two-layer and multi- 
layer models have been discussed by Wait [1958, 
1962]; Madden [1961]; Jean et al. [1961], and a par- 
ticularly successful two-layer model has been con- 
structed by Chapman and Jones [1964a, b]. This 
model, which is indicated in figure 5, has been applied 
in the frequency range of earth-to-ionosphere cavity 
resonances and also for frequencies up to 1 kc/s. 

Such models are usually derived by trial and error 
procedures and they are difficult to correlate with con- 
tinuous ionosphere profiles, such as shown in figure 4. 

3.5. Models of Continuously Varying Ionosphere 

The isotropic exponential model gives a reasonable 
approximation to the ionospheric conductivity versus 
height profile through the lower ionosphere layers, and 
the propagation parameters of the model can be deter- 
mined from closed form expressions of the surface 
impedance [Galejs, 1961a, 1962a, b, 1964a]. Calcu- 
lated resonance frequencies and Q factors that are 
based on (3), (4), and available values of S [Galejs, 
1962a] are listed in tables I and II. Calculations are 
based on daytime or day and nighttime ionosphere 
models. The Q factors considered also the effects 
of cosmic ray ionization [Bourdeau et al., 1959] that 
give a nearly exponential conductivity versus height 
variation as indicated by Galejs [1962a]. Cosmic ray 
ionization decreases the effective Q figures, in particu- 
lar for the lower frequencies. The calculated Q values 
are still somewhat higher than those estimated from 
experimental data. 

TABLE 1 
RESONANCE FREQUENCIES OF THE SPHERICAL SHELL BETWEEN THE EARTH AND IONOSPHERE 



■ 




12 3 4 


1 ^ c 

O z 


LOSSLESS CAVITY 




10.6 18.3 25.9 33.5 


DAYTIME WITH h = 45 KM 




7.6 13.5 19.4 25.4 


AVERAGE DATA OF DAY WITH h = 45 KM 
AND OF NIGHT WITH h = 60 KM 




8.05 14.25 20.5 26.8 


MEASUREMENTS 




8 14.1 20.3 26.4 



Q FACTORS OF EARTH- IONOSPHERE CAVITY RESONANCES 



f IN C/S 


10 30 100 


o 

H 

u 


DAY WITH h = 45 KM, NIGHT WITH h = 60 KM 


8.5 7 6.7 


DAYTIME ONLY WITH h = 45 KM 


6.5 5.4 5.2 


DAY WITH h = 45 KM, NIGHT WITH h = 60 KM, 
COSMIC RAY IONIZATION 


5.5 6 6.4 


MEASUREMENTS 


4 6 



1046 



3.6. Magnetic Field Effects 

The magnetic field effects for a homogeneous sharply 
bounded ionosphere model have been discussed by 
Wait [1962] in a quasi-longitudinal approximation for 
a superimposed radial magnetic field and also for a 
transverse magnetic field (east-west or west-east propa- 
gation along the magnetic equator). 

The propagation parameters have been determined 
by Galejs and Row [1964] and Galejs [1964a, 1964b] 
in the presence of a horizontal magnetic field. These 
data provide an indication of propagation character- 
istics of the waves near the magnetic equator, but 
they do not apply directly to the earth-ionosphere 
resonance problem, where the H vector is nearly hori- 
zontal only over a small surface portion of the earth. 

The earth-ionosphere cavity resonance frequencies 
have been considered in the thesis of Thompson 
[1963] for a radial magnetic field. He uses a multi- 
layer ionosphere model that approximates the con- 
tinuous altitude variation of the conductivity tensor. 
The azimuthal field variation is expressed in terms of 
the Legendre polynomials P n (cos 6) for the nth mode. 
The wave equation has four different solutions for 
each of the ionosphere layers, but it has only two 
solutions for the space between the ionosphere and 
the perfectly conducting ground surface. The solu- 
tion of the problem and the resonance condition for 
the mode n are obtained after multiplying a sequence 
of 4 X 4 matrices. In these computations the displace- 
ment currents of the ionosphere are neglected. The 
computed resonance frequencies are in agreement 
with measurements, but the Q values are in the range 
of 15 to 20, with the exception of nighttime anisotropic 
solutions, where Q = l to 10. The lowering of the 
nighttime Q figures by the magnetic field has been 
attributed to wave penetration through the ionosphere 
and to energy thus escaping from the earth-ionosphere 
cavity. 

Another solution to the problem of earth's 
ionosphere cavity resonance in the presence of 
a radial magnetic field is given by Galejs [1965]. 
The azimuthal field variation is expressed in terms of 
the Legendre functions P v (— cos 0), and the solution 
of the wave equation is set up so that the four solutions 
of the ionospheric layers where a > w€ degenerate 
into two solutions for the lower ionospheric layers 
where a < < co6q and for free space. After comput- 
ing the surface impedance at the lower boundary of 
the ionosphere, the propagation parameters are deter- 
mined from the solutions of the modal equation for TM 
modes, althought the coupling between TE and TM 
modes is considered in the processes of computing the 
surface impedance. Only the vertically polarized TM 
modes are assumed to propagate below the ionosphere. 
Such assumptions have been justified for a single layer 
ionosphere model at higher frequencies [Wait 1962, p. 
269] and for a thin ionospheric shell also for fre- 
quencies in the resonance range [Wait, 1964c]. 

In the numerical calculations the anisotropic 
ionosphere models are those of figure 4 and further 
calculations are made for ionosphere models which 



consider electrons and ions at an operating fre- 
quency of 20 c/s, as indicated in figures 3 and 4 
of Galejs [1964b]. In the isotropic ionosphere 
models cr ( ) remains unaltered, en = <x , and o-> —* 0. 
The calculated resonance frequencies f n are listed in 
table III. The isotropic day or night models exhibit 
too high resonance frequencies, but the average for 
the anisotropic day or night models gives nearly 
correct results, although the first resonance occurs 
near 7.7 c/s. At night the various anisotropic models 
give different results, while their results are nearly the 
same for daytime or for the isotropic models. The 
calculated Q factors are listed in table IV. At daytime 
the Q factors are lowered slightly due to anisotropy. 
The losses are principally due to absorption in the 
lower D region where the collision frequencies are 
high and the magnetic field has only small effects. At 
nighttime the anisotropy reduces the Q values dras- 
tically, and the ions have a considerable effect. For 
the isotropic models the energy remains below the 
ionosphere, the losses in the lower ionosphere are small 
and the Q values are high. For the anisotropic models 
some energy penetrates the ionosphere and escapes 
from the earth to ionosphere cavity. Hence the Q 
figure is low. 
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The previous calculations have assumed a constant 
radial magnetic field (S pole) all over the globe. A 
reversal of the static magnetic field (N pole) causes 
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negligible differences in the propagation parameters or 
in the values of f n and Q. Hence a model using the 
average of the radial magnetic field may provide a 
first order estimate of the anisotropy effects. The 
present calculations assume the radial magnetic of the 
polar regions to be extended over the whole surface of 
the globe and the anisotropy effects should be too 
pronounced. The analytical results are highly 
dependent on the detailed structure of the nighttime 
anisotropic ionosphere, and further studies of it may 
be in order. 

4. Source Distributions and Noise Spectra 

Most of the ELF energy is of terrestrial origin and is 
caused by thunderstorm activity, although there are 
some high latitude events of an apparent extra- 
terrestrial origin. Such emissions of approximately 
800 c/s have been detected near Kiruna, Sweden 
[Egeland, 1964] and of 14 to 17 c/s on Kerguelen Island 
in the Southern part of the Indian Ocean [Gendrin and 
Stephant, 1964, private communication]. The 
Spectran recording of the signals received at Kerguelen 
is shown in figure 6. The upper trace shows a normal 
signal with resonances near 8, 14, 20, and 25 c/s. 
There is an additional signal of 14 to 17 c/s on the two 
lower traces. When neglecting the presence of such 
extraneous signals the terrestrial noise spectra can 
be calculated after estimating the dipole moments of 
the sources and using propagation parameters that are 
appropriate to the selected ionosphere model. 
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Figure 7. Estimated lighting flash spectrum. 



Raemer [1961a, b] has estimated the power spectrum 
of the dipole moment in a median lightning flash Ids (or 
the first derivative of the electrostatic moment dMjdt) 
based on statistics of lightning strokes compiled by 
Williams [1959]. Later work of Raemer [1963 private 
communication] has lead to a revised estimate. These 
spectra are shown in figure 7. Pierce [1963] has com- 
puted the magnitude of the frequency spectra of the 
second derivative d 2 M/dt 2 , and Rycroft [1963] has 
calculated the frequency spectra of dMjdt using a 
similar representation of the lightning moments. The 
square of the frequency spectrum computed by Rycroft 
[1963] is also indicated in figure 7. 

An equivalent power spectrum for the dipole moment 
of terrestrial noise sources has been deduced by 
Harris and Tanner [1962] from the measurements of 
Balser and Wagner [1960] by estimating the complex 
propagation constant within this frequency range in 
a "trial and error" approach. Harris and Tanner 
[1962] do not require an ionosphere model or any 
knowledge of lightning waveforms and their statistics 
for obtaining this equivalent power spectrum. 

After determining the source spectra and the propa- 
gation parameters the power spectrum of the received 
noise may be computed from 



G(ioj) 



g{i<o)\v(p+l)\ 



8(ak h) 2 | v + \ | [cosh (2tt Im v) - cos (2tt Re v) ] 
sinh[2Im^(7r-fl)] 
(- Im v) 

sin[(2Re^ + l)(7r-fl)-(7T/2)] ]fl 2 



Re^+ \ 



(9) 



FIGURE 6. Spectran recordings of noise waveforms at Kerguelen 
Island. 



where the v is given by (2) and where the sources are 
assumed to be uniformly distributed in the 6 interval 
of S\< < 02- The derivation of (9) starts from (1) 
and Legendre functions P v {~ cos 6) are replaced by 
an asymptotic expansion. There are also equivalent 
spectral representations in terms of Legendre poly- 
nomials [Raemer, 1961a, b]. 



1048 



4.0 

CO 

I 3 - 2 

> 

5 2.4 



g'- 6 



1.8 " 



01 L_ 






THEORETICAL 

EXPERIMENTAL 



Figure 8. 



16 20 24 

FREQUENCY- c/s 



36 



ELF Noise spectrum— effective ionospheric height 
assumed constant. 
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The noise spectrum of Raemer [1961b] has been 
used for computing the observed noise spectra in 
conjunction with a sharply bounded ionosphere model. 
This spectrum is shown in figure 8 and is seen to exhibit 
excessively damped higher resonance peaks. 

The exponential ionosphere model [Galejs, 1961b, 
1962a, 1964a] produces a better agreement with meas- 
urements as seen in figures 9 and 10. In this work 
propagation is assumed to be uniform around the sur- 
face of the earth and the propagation parameters 
represent the average of day and night conditions. 
The sources are assumed to be uniformly distributed 
over the 6 interval indicated in the figures (0 is denned 
with respect to the observer). 

5. Diurnal and Seasonal Changes of 
Power Spectra 

Balser and Wagner [1962a] have recorded diurnal 
variations of measured noise spectra. Their meas- 
urements of the diurnal variations of the power level 




UNIVERSAL TIME 



Figure 11. Diurnal variation of the resonance modes at Ipswich, 
Mass., February 1961 (vertical electric field), Chambon-la-Foret 
(near Paris) July 1962 and at Tromso (Norway) April 1962 (hori- 
zontal magnetic field). 



at the first three resonance frequencies are shown in 
figure 11 together with similar measurements by 
Gendrin and Stefant [1962a]. Diurnal variations of 
the fields near the two lower resonance frequencies 
have been measured at Kingston, R.I., by Polk, Huck, 
and Yu [1963, private communication]. The maxi- 
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mum signal intensity occurs between 1500 and 1800 
UT. In an amplitude versus time plot the signal peak 
is broader during summer than winter. Seasonal and 
diurnal variations of the noise power in the ELF fre- 
quency band below 35 c/s have been described by 
Wright [1963] without singling out the individual 
resonance modes. The measurements are made at 
Byrd Station, Antarctica (80°S, 120°W) and show peak 
intensities near 1000 UT during spring and summer 
and near 2300 UT during most of the year. 

The gross features of the spectrum, such as the 
number of peaks and their approximate time of occur- 
ance, may be explained qualitatively by noting areas 
of globe with particularly high thunderstorm activity 
and by measuring the polar angle from these source 
regions to the observation point. The vertical elec- 
tric field E r and the horizontal magnetic field H<t> of 
the mode n exhibit the dependence of 



E r ~Pn(cOs0) 

H^-~P n (cos6\ 



(10) 
(11) 



where P n (cos 0) is a Legendre polynomial of the order 
n. A particular geographical region will contribute 
to the measured values of E r or H& in the mode n if 
the right hand sides of (10) or (11) are sizeable for those 
values of 6. Such considerations have been reported 
by Balser and Wagner [1962a] and Gendrin and 
Stefant [1962a]. 

More detailed calculations have been carried out by 
Abraham [1964, private communication; Galejs 1964a] 
who used available estimates of geographical distribu- 
tions of the noise sources and who calculated the day 
and night propagation effects after expressing the fields 
as a sum of recirculating azimuthal waves. 

In the expression for the electric field component 
E r of (1), the Legendre function is replaced by the 
leading term of its asymptotic expansion. The func- 
tion sin vir is expanded into 



(sin v7T)~ 1 =-2i j? e« 2n + 1 



)VTT 



(12) 



which is valid for Im v > 0. This results in [Abra- 
ham, 1964, private communication] 



E r - 



Ids 



v(v+ 1) 



|3f 



2(oeha 2 V2tt(i/+ J)sin0 



i(v + \)(0 + 2nir)-i2n- 



Hv+ l)[2ir(n + l)-0}-i{2n+l) 



(13) 

The first term of the summation represents a wave 
which has traveled the direct distance from the source 
to the observation point (0) plus n full circles around 
the sphere (2m7t), while experiencing a phase shift of 
(— 7i/2) radians at each traversal of the source point 
or its antipode. The other term of the summation 
represents a wave which has traveled toward the ob- 
servation point across the antipode of the source 
(27T — 0) and has made furthermore n circles around the 



sphere {2mr) and has also experienced a phase shift 
of (— 7r/2) radians at each traversal of the source point 
or its antipode. This physical interpretation of (13), 
which is equivalent to (1) for v = const., permits a 
heuristic consideration of the v differences between the 
day and night hemispheres. The multiplicative v fac- 
tors in front of the n summation are assumed to be con- 
stant and equal to v which is intermediate between 
the day and night v values. The 6 factors of the ex- 
ponentials are replaced by integrals in order to 
account for the v variation seen by the circular wave 
front as it propagates in the ^-direction with respect 
to the source. The amplitude changes at the day 
and night boundary can be also considered for each 
traveling wave by assuming that the energy of the 
wave remains constant after passing the day and night 
boundary. This results in decreased amplitudes for 
the larger ionospheric height of the night hemisphere. 
In the numerical work the surface of the earth was 
subdivided into 10° cells. Each cell had a constant 
surface density of noise sources but a variation was 
allowed from cell to cell. The plots of thunderstorm 
days over the world were used to set relative values 
of lightning activity for each of the cells, by assuming 
that the distribution of squared dipole moment per 
unit area is directly proportional to the thunderstorm 
day plots according to the Handbook of Geophysics 
[I960]. The diurnal variations were accounted for 
by a multiplicative factor in terms of local time at 
the source [Williams, 1959]. The final equation 
used in the spectrum calculations is similar to (9), 
although it is algebraically more involved [Abraham, 
1964, private communication]. 

The diurnal variations of the noise power at the 
first three Schumann resonances have been computed 
using the exponential ionosphere model of figure 5 
and are shown in figure 12. The calculations have dis- 
regarded the noise sources near the antipode where 
the usual asymptotic approximation of the Legendre 
functions is inaccurate. Figure 12 shows a qualitative 
agreement with the measurements of Balser and 
Wagner [1962a], that are shown in figure 11. 

An example of seasonal power variations is shown 
in figure 13. It indicates a decreased noise power 
near the higher frequencies and/or during winter 
months. 

Based on measurements at one location of the earth 
it is also possible to predict the noise levels at other 
locations using this procedure of noise spectrum 
calculations. The noise that is due to a worldwide 
distribution of thunderstorms depends on the fre- 
quency variation and also on the absolute value of 
the power spectrum of the mean dipole moment of 
the source g(i(o, 0, (/>). The and (j) dependence may 
be estimated from the thunderstorm frequency maps 
as was indicated before, which leaves the absolute 
value and the frequency dependence of giico) as an 
unknown. Assuming that g(ico) is the same for all 
noise sources, it may be determined by making the 
calculations to agree with measurements at one par- 
ticular location. This gives a uniquely determined 
noise distribution over the surface of the earth and 
the accuracy of this model can be checked by com- 
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FIGURE 12. Calculated diurnal variation of noise power at the 
peak of the resonance modes — wintertime — Boston, Mass. 
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FIGURE 13. Calculated noise spectrum for local noon in northern 
hemisphere. 



FIGURE 14. Comparison of calculated field strength's data with 
July 1963 measurements by DECO and RCA. 

paring the calculations with measurements at other 
locations. In an example, measurements by DECO 
at Boulder, Colo., [Maxwell, 1963, private communi- 
cation] and by RCA at Long Island, N.Y. [Powers, 
1964, private communication], were used to determine 
g(ico). The calculations by Abraham [1964, private 
communication] based on these values of g(i(o) have 
been compared with simultaneous measurements 
made in Malta and Alaska [Maxwell, 1963, private 
communication]. As seen from figure 14 the analytic 
curves agree in shape with the measurements, but the 
normalization with respect to Boulder gives somewhat 
higher predictions of the noise level. This may pos- 
sibly be caused by summertime thunderstorm activity 
in the mountains near Boulder, Colo., that would give 
significant near-fields and a high estimate for the func- 
tion giico). The calculations are extended above the 
region of the earth-ionosphere cavity resonances to 
show that more significant noise level changes occur at 
higher frequencies. 

In the calculations leading to figures 12 through 14 
the day and night effects are considered most ac- 
curately if the day and night boundary is located sym- 
metrically with respect to source (i.e., source at local 
noon or midnight). For nonsymmetrical day and 
night boundaries it is difficult to estimate the effects 
of wave front distortion. When neglecting such distor- 
tions (or assuming propagation strictly in the 6- 
direction), this method can account for v{6) variations, 
which may be caused by ionosphere parameter changes 
or by changes of the static magnetic field along the 
surface of the earth. The assymptotic expansions of 
P v (—cosd) used in (13) or (9) can be modified to 
increase their accuracy for observations near the 
source or its antipode [Abraham, 1965]. 

6. Variations of Resonance Frequencies 

Diurnal frequency variations of the earth-ionosphere 
cavity modes have been reported by Balser and Wagner 
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[1962b], Gendrin and Stefant [1963], Rycroft [1963] 
and by Chapman and Jones [1964b]. The resonance 
frequencies vary by about ±0.2 to 0.3 c/s, but the 
changes are not simultaneous for the different modes; 
a low resonance frequency of one mode occurs simul- 
taneously with a high frequency of a different mode. 
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Figure 15. Diurnal frequency variations measured by Gendrin 
and Stefant. n 
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Figure 16. Diurnal frequency variations measured by Rycroft. 
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FIGURE 17. Nuclear effects on earth-ionosphere cavity resonances. 

Measurements by Gendrin and Stefant. 



A high altitude nuclear explosion has been shown to 
affect the earth-to-ionosphere cavity resonances by 
simultaneously lowering all the resonance frequencies 
[Gendrin and Stefant, 1962a and b; Balser and Wagner, 
1963]. Tepley et al., [1963] noted also an increase of 
signal strength for observations near the conjugate 
point of the explosion. The lowering of resonance 
frequency shown in figure 17 has been attributed to 
worldwide lowering of the effective ionosphere in a 
model of a single layer isotropic ionosphere [Gendrin 
and Stefant, 1962b]. The depression of the resonance 
frequency can be also explained by assuming an ex- 
ponential daytime ionosphere to be effective around 
the globe [Row, 1964, private communication]. As an 
illustration the daytime resonance frequencies of 
table 1 compare with the measured resonance fre- 
quencies of figure 17. There are VLF observations 
[Zmuda et al, 1963] that can be explained by a lower- 
ing (or increasing electron density) of the Z)-region 
boundary. Brady et al., [1964] interpreted VLF di- 
urnal phase shifts to show that some enhanced D- 
region ionization persists for several days following 
the event. 

The lowering of the resonance frequency following 
a nuclear explosion has been attributed by Chapman 
and Jones [1964a] in their two-layer ionosphere model 
solely to a decreasing conductivity of the £-region. 

After measuring a change of the resonance frequen- 
cies following a nuclear explosion, it is possible to 
advance simple ionosphere models for explaining a 
particular change. However, no attempt has been 
reported to predict resonance shifts for different 
nuclear yields and detonation altitudes. This rather 
complex problem requires a quantitative evaluation 
of the various nuclear mechanisms that may lead to 
ionosphere perturbations. If such information is 
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available, the computation of the resonance shifts 
is relatively simple. It would be also of interest to 
investigate the effects of sudden ionosphere dis- 
turbances (SID) and to discriminate between SID's 
and nuclear effects. 

The resonant frequency changes can be considered 
more accurately with nonuniform models of the propa- 
gation geometry. A perturbation method [Wait, 
1964b] or numerical techniques [Madden and Thomp- 
son, 1964] have been used for solving such problems. 

Madden and Thompson [1964] compute the local 
surface impedance and propagation constant for a 
uniform cylindrical shell between the anisotropic 
ionosphere and the earth using matrix multiplication 
techniques. There is a smooth variation of the 
parameters between polar and equatorial regions 
and also between day and night hemispheres. These 
local surface impedances and propagation constants 
are assumed to characterize a spherical transmission 
line of nonuniform parameters. The resonance fre- 
quencies and Q factors are computed numerically 
using a lumped parameter approximation of the two- 
dimensional transmission line. This model has been 
applied to explain the diurnal frequency variations 
and also the frequency changes due to a nuclear 
explosion. 



7. Areas of Future Work 

There are numerous observations of earth-iono- 
sphere cavity resonances and existing ionosphere and 
source models account for most observed charac- 
teristics like the resonance frequency, Q factors, 
shape of the noise spectra, and diurnal variations of 
the received power. 

Future work appears to be required in a systematic 
study of diurnal and seasonal amplitude and frequency 
variations of resonances. These changes should also 
be correlated with ionospheric disturbances. The 
parameters of the nighttime ionosphere require also 
a future clarification. The variations of the ionosphere 
parameters and the variations of the magnetic field 
vector along the surface of the earth should be in- 
corporated in future calculations [Wait, 1964b; Mad- 
den and Thompson, 1964]. Much work remains to be 
done for assessing nuclear effects on the earth iono- 
sphere resonance and for discriminating between 
lower level nuclear and sudden ionospheric disturb- 
ances. 



lossy walls is 
J_ = 

Qi 



J J E x H* • ds 



tfSSSE 



E^+^J / /H • U*dv 



(14) 



where the surface integral is extended over the cavity 
walls and the volume integral covers the inside of the 
cavity. The ground surface is assumed to be a perfect 
conductor and Ee is related to H^ at the ionospheric 
boundary by the surface impedance Z s as Eq = Z s H^. 
However Ee is assumed to be small relative to E r for 
vertically polarized waves. Equation (14) is rewritten 



1 = 2 f J ReZslff^q + A) 2 sin dddd± 
Qi (o / / / [e |£r| 2 + /xo|//<>| 2 ]r 2 dr sin 6d6d(j> 



2ReZ, 



J \H*\ 2 sin Odd 
(oh J[€o|£r| 2 + jU, |H4| 2 ] sin 9d0 



(15) 



(16) 



where a is the radius of earth, h is the ionospheric 
height and where the cavity is assumed to be uniform. 
For excitation by a radial electric dipole of moment 
Idr the field component E r and H^ follow from (21), 
(61), and (62) of Galejs [1964a] as 



E r = 



ildr n(n + l)(2n+l) 



4>7Tho)€a 2 



D n 



P„(cos 0) (17) 



H^ — — 



Idr (2/1+1) d 



kirha 



D n Te p » {cose) 



(18) 



with D n = n{n+l)-v{v+\). Substituting (17) and (18) 
in (16) and evaluating the integrals gives 



Qi 



2ReZ, 



(O/JLoh 



n(n+l) 



(koa) 2 



+ 1 



(19) 



where ko=coy/JLoeo . From the approximate modal 
equation 



'=yjl + 



iZ s 

(O/JLoh 



it follows that 



2ReSImS = 



ReZ s 

(O/JLoh 



(20) 



(21) 



The work of Sylvania reported in this paper has Applying (21) and the resonance condition 
been supported principally by the Office of Naval 
Research under Contract Nonr 3185(00). n(n+ 1) ~ (k n a ReS) 2 



(19) simplifies to 



8. Appendix. Q Factors 

The basic definition of the Q factor of a cavity with 



1 4ReSImS 
Qi (ReS) 2 + l 



(22) 



(23) 
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which was listed as (4). When considering the dif- 
ferences of E r , //(£, Z s , and h between the day and 
night hemispheres it follows that 



1 2[ReS ImS| day + ReS ImS| night ] 



Qi 



(ReS 



! +l 



(24) 



The form (5) or (6) is obtained by neglecting the E 
integral or H integral, and by doubling the value of 
the H integral or E integral in (14). 
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